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We study dilute aqueous solutions of charged disk-like mineral particles (Laponite) by a combina- 
tion of static light scattering (SLS) and small-angle x-ray scattering (SAXS). Laponite solutions are 
known to form gels above a certain critical concentration that must be described as non-equilibrium 
states. Here we focus on the investigation by SLS and SAXS at concentrations below gelation 
(c < 0.016 g/L) and at low concentrations of added salt (0.001 and 0.005 mM). Thus, we have ob- 
tained the scattering function of single Laponite platelets as well as the structure factor describing 
their interaction at finite concentration. A detailed analysis of the combined sets of data proves 
that the solutions are in a well-defined equilibrium state. Moreover, this analysis demonstrates the 
internal consistency and accuracy of the scattering functions obtained at finite concentrations. We 
find that Laponite particles interact through an effective pair potential that is attractive on short 
range but repulsive on longer range. This finding demonstrates that Laponite solutions exhibit only 
a limited stability at the concentration of added salt used herein. Raising the ionic strength to 
0.005 mM already leads to slow flocculation as is evidenced from the enhanced scattering inten- 
sity at smallest scattering angles. All data strongly suggest that the gelation occurring at higher 
concentration is related to aggregation. 

PACS numbers: 61.10.-i, 61.25.-f, 64.70.Ja 
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I. INTRODUCTION 



Laponite consists of disk-like clay particles with a 
thickness L « 1 nni and a radius R w 12.5 nm. Each 
platelet carries a few hundred elementary charges. For 
a long time, aqueous suspensions of Laponite have been 
extensively investigated as a model system for disk-like 
colloids Q, 1^ 0|- A theoretical description of such a 
well-characterized model clay appeared feasible. Early 
attempts focussed on the structure of the electric dou- 
ble layer around a single platelet 0, IE 0) or around 
two parallel platelets Much of the recent work was 
initiated by the theoretical treatment by Hansen and 
coworkers [3- In the model of Hansen et al. Laponite 
platelets and their associated electric double layers are 
represented by nonintersecting discs carrying a constant 
electrostatic quadrupole moment. The spatial arrange- 
ment of these objects has been investigated by Monte 
Carlo simulations. The results of these studies point to 
a reversible sol-gel transition at platelet number density 
n > 0.06 (corresponding to a platelet concentration 
of 3.4 % by weight), where the structure of the gel phase 
is strongly reminiscent of a "house-of-cards" jlfll- How- 
ever, the quadrupolar disc model has turned out to be 
an oversimplification. In particular it is inadequate at 
very short range, where any multipolar expansion breaks 
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down, and at long range, where electrostatic interactions 
between the platelets are exponentially screened. 

A screened electrostatic potential between two arbi- 
trary oriented, charged platelets has been worked out 
within the framework of a linearized Poisson-Boltzmann 
theory P, [ij, [ij, III . At a fixed center-to-center dis- 
tance, this electrostatic potential (always repulsive) is 
maximized for co-planar platelets, which corresponds to 
the maximum overlap of electric double layers, and min- 
imized when the platelets are co-axial and parallel. An 
intermediate potential is found for T-shape perpendicular 
platelets. At present only preliminary results of Monte 
Carlo simulations for platelets interacting through such 
an anisotropic potential are available, indicating that the 
platelet fluid is in the isotropic phase for platelet number 
density n < 0.1 R'^ fl^. 

Later an interaction site model of Laponite was put 
forward, and investigated by molecular-dynamics simula- 
tions [l3| and integral equations ^IS] . In this model each 
platelet carries discrete charged sites rigidly arranged on 
an array inscribed in the disc. Sites on different platelets 
interact via a repulsive screened Coulomb interaction. 
Upon varying the platelet number density and screening 
length, a sol, gel, and crystal phase have been identified. 
Additional rim charges of opposite sign have been found 
to lead to T-shaped pair configurations and clustering of 
the platelets. It should be noted, however, that these 
models treat the Laponite particles as uniform. 

Summing up, several mechanisms have been proposed 
as origin for the gel formation in Laponite solutions by 
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now: First of all, the classical answer to this problem is 
the formation of a "house of cards" -structure by attrac- 
tive interactions between the platelets, possibly caused 
by positive charges on the rims 10] . This includes more 
recent approaches taking into account the charge distri- 
bution around the platelets [T^. 

As an alternative to this modeling, it was argued that 
this type of aggregation will only take place at higher 
ionic strength (> 0.01 mM) and may be ruled out as 
cause for gelation. Instead, the repulsive interaction 
between the c harg ed platelets leads to a glassy state 
H m 111 in Efl m pi m m m m. Slow time- 
dependent processes within the gel were interpreted in 
terms of an aging of the glass. Very recently, it has been 
hypothesized that Laponite suspensions behave like an 
attractive glass ||23, |23 . A comprehensive discussion of 
the treatment of possible glassy states in Laponite sus- 
pensions including a comparison with experimental data 
was recently given by Mongondry et al. |29j. 

A markedly different model recently proposed by Nico- 
lai et al. invokes aggregation in order to explain the pro- 
cess of gelation ,31j . Weak attraction may lead 
to a slow aggregation of the particles even at low ionic 
strength. As a consequence of this, no well-defined phases 
or phase transitions can be identified anymore. The ob- 
servations made by Nicolai et al. |2^ are in agreement 
with a study of Martin et al. |3^. Hence, the origins 
of the sol-gel transition have not yet been fully clarified 
despite of more than two decades of research. 

Most of the theoretical expositions on gelation and 
structure formation in Laponite suspensions have as- 
sumed a system of monodisperse platelets, i.e., thin disks 
with a radius of 12.5 nm and a thickness of 1 nm. How- 
ever, it has become possible to visualize directly the 
Laponite particles by various techniques as AFM or elec- 
tron microscopy. From these micrographs it is obvi- 
ous that La pon ite exhibits a considerable polydisper- 
sity 0, I35I Isd W^ . A second point that is common 
to many experimental study is the analysis of obvious 
non-equilibrium states reached after passing the gel line. 
As the matter of fact, considerable efforts have been de- 
voted to the investigation of the exact position of the gel 
line. The experimental difficulties of the investigations 
of such time-dependent states are at hand. Moreover, no 
clear information on the nature of the interparticle po- 
tential can be derived from essentially non-equilibrium 
states. An advanced understanding of Laponite suspen- 
sions, however, requires a clear delineation of the range 
of repulsion together with a clear assessment of a possible 
attractive part of the potential that may come into play 
at short distances [Tol| . 

In this paper we present a comprehensive study of 
Laponite suspensions in the dilute regime by static 
light scattering (SLS) and small-angle x-ray scattering 
(SAXS). We purposely avoid higher concentrations lead- 
ing to gelled states. Instead, particular care is devoted to 
the proof that we do look at equilibrium states indeed. 
Hence, data are taken at different small concentrations 



of the solute and extrapolated to infinite dilution. This 
leads to the scattering intensity Iq (q) of single platelets 
as the function of q, the magnitude of the scattering vec- 
tor q {q — (Att/X) sin(6'/2); A: wavelength of radiation; 9: 
scattering angle). Ia{q) in turn can be modelled to give 
the polydispersity. The scattering intensity I{q) obtained 
at finite but small concentrations is used to determine 
the "measured structure factor" Spiiq) which is related 
to the pair correlation function g{r) and the interaction 
potential V{r). However, the appreciable polydispersity 
of the Laponite particles precludes a direct analysis of 
Sniq) to yield V{r). Here we shall address this problem 
within the framework of the " Polymer reference interac- 
tion site model" (PRISM) . Using the approach developed 
previously for suspensions of bidisperse platelets Is^ , size 
polydispersity of the platelets will be taken into account 
by a multicomponent PRISM integral equation theory. 
We shall demonstrate that the combination of precise 
scattering data with the multicomponent PRISM theory 
leads to an unambiguous assessment of the interaction 
potential in solution. 

The paper is organized as follows: After the section 
Experimental we shall delineate in section HI the theory 
of light scattering of dilute suspension. In particular, the 
check of the measured data for consistency is given. This 
section contains also the description of the multicompo- 
nent interaction site model used for calculating structure 
factors. Section IV will give describe the comparison with 
measured scattering data and the resulting effective in- 
teraction potential. A final brief section V will conclude 
this paper. 



II. EXPERIMENTAL 
A. Materials 

Laponite RD was obtained from Laporte Industries 
Ltd. and is used for our study without further purifica- 
tion. Laponite powder was added to a 10^"^ mol/1 NaCl 
solution in distilled water and mixed at high speed with 
a magnet stirrer for 4 hours. The dispersions were ini- 
tially turbid, but became clear in 15-20 minutes, indicat- 
ing that the particles were fully hydrated. All Laponite 
samples (volume fraction ranging from 0.02 % to 0.16 %) 
were prepared separately, without using a stock solution. 
Vigorous stirring of the solutions was found to be nec- 
essary for dispersing particles homogeneously. SAXS- 
experimcnts demonstrated that incomplete dissolution 
leads to the formation of clusters and large aggregates. 

The pH values of all the samples were adjusted to 10 
by using NaOH in order to avoid degradation or disso- 
lution of Laponite particles. Dispersions were kept still 
overnight before any use in order to achieve equilibrated 
state. For SAXS experiments the samples were measured 
without further treatment, as no difference was observed 
between the scattering intensities of unfiltered dispersion 
and filtered one by using 0.22 mm filters (Millipore). For 
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static light scattering experiments the dispersions were 
filtered by using 0.22 mm filters (Millipore) into quartz 
cuvettes to eliminate dust and remaining aggregates. 



B. Methods 

The scattering intensity I{q) was obtained by static 
light scattering in the region of smallest g-values and by 
SAXS for higher g-values. Static light scattering were 
performed by a Sofica SLS-Spectrophotometer, using ver- 
tically polarized light with a wavelength of 632.8 nm. 
The sample is contained in a cylindrical 10 mm diameter 
quartz cuvette, which is immersed in an index-matching 
toluene bath of 80 mm diameter. Scattered light is de- 
tected by a photomultiplier, which is mounted on a go- 
niometer arm. The angle range is 35-125 deg., corre- 
sponding to q range of 7.9 x 10""^ — 2.3 x 10^^ nm^^. 

For SAXS experiments, the data are composed by 
two sets in different q ranges. One set of data in rel- 
atively lower q range were measured at the European 
Synchrotron Radiation Facility (ESRF) in Grenoble on 
the high brilliance beam line ID2. The instrumental set- 
up for Laponite samples is as following: wavelength A = 
0.0995 nm with energy 12460 eV, sample-detector dis- 
tance d = 10 m which covers an experimental range 
in q between 2.15 x 10"^ to 0.67 x 10"^ nm'^ A 
two-dimensional image-intensified charge-coupled device 
(CCD) is employed as the detector, which allows the 
study of the (an)isotropy of the scattering. The other 
set of SAXS data in q range from 0.07 to 6 nm~^ were 
obtained in our lab by a Kratky commercial camera. The 
cameras and the data treatment have been described else- 
where 39] . 

Density measurements were done in dilute solutions us- 
ing a DMA 60 apparatus supplied by Paar, Graz, Austria. 
This gave a mass density of 2.42 ± 0.12 g/cm'^ for the dis- 
solved Laponite particles. The refractive index increment 
was measured using the Refraktometer alpha-Ref Typ 
1 apparatus supplied by SLS-Systemtechnik (Freiburg, 
Germany) operating at the wavelength of the apparatus 
using for static light scattering (632.8 nm). 



The scattering intensity of a single platelet Iq (q) may be 
used to define the form factor P{q) that is normalized to 
unity at q = by 



Io{q)^Vp\p-p^)'P{q), 



(2) 



where Vp denotes the volume of the particles, j? is the 
scattering length density of the particles and pm is the 
respective quantity of the solvent (see the discussion of 
this point in Ref. [S Q). Since the volume fraction 
is given by ^ = {N/V)Vp it follows that I{q = 0)/(j) = 
Vpip — Pm)^ where p — pm is the contrast of the particles 
in the respective solvent. If the concentrations of the 
solute is expressed through the weight concentrations c, 
Vp must be replaced by the molecular weight A/, that is 



liq) = KcMP{q)S{q) . 



(3) 



In case of SAXS the constant K then follows from the 
number of scattering units per unit mass of the solute. 
For light scattering, K is given by 



4^2^)2 



K = 



(4) 



where no is the refractive index of the solvent, (dh/dc) is 
the refractive index increment, A denotes the wavelength 
of the light, and Na is Avogadro's number. 

With comple te g enerality, P{q) may be expanded for 
small q to give |43, 



P{q) « exp(- 



-Rgq 



(5) 



where Rg is the radius of gyration of the particles. Io{q) 
may hence be extrapolated to q — 0. In general, a platelet 
can be modelled as a circular disc of radius R and thick- 
ness L. For randomly oriented monodisperse platelets we 
have lis 
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P{q) = j sin?9 





4 sin((ji /2 cos d) Ji [qR sin d) 
q^RL cos I? sin I? 



(6) 

where Ji denotes the cylindrical Bessel function of first 
order. 



III. THEORY 

A. Analysis of particle interaction by light and 
small-angle x-ray scattering 

We consider a dilute isotropic suspension of N 
monodisperse platelets per volume V . The scattering 
intensity I{q) can be rendered as the product of the in- 
tensity I[j{q) of single platelets and the structure factor 
S{q) through [Ulllll^ 



/(g) = ^/o(g)5(g). 



(1) 



B. Interaction in dilute solution 

The dependence of S{q) on concentration may be ex- 
pressed in the usual virial series by flO, 41« |4^ 



S{q) 



1 + 2Bi^pp(f) + 3Capp(i 



(7) 



where the apparent virial coefficient -Bapp depends on q 
through [iSlillil 



^app — 4 



(^/6)d; 



eff 



(1 



1 j2 2 



(8) 
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and def / is a measure for an effective diameter of inter- 
action in the dilute regime. The quantity Capp is the 
correction related to the third virial coefficient. Eq. © 
can be derived easily from the fact that binary interac- 
tion prevails if the concentration is low enough. There- 
fore the pair correlation function g(r) follows simply by 
g{r) = exp[—V{r)/kT] ji^.The scattering data obtained 
in the dilute regime can therefore be extrapolated to van- 
ishing concentration by plots of (j)/I{q) against concen- 
tration. Moreover, Eq. © suggests that the slope of 
these plots that give the apparent virial coefficient Bapp 
will diminish with increasing value of q |47j . Hence, a 
precise analysis of the scattering intensity I{q) obtained 
for different concentrations leads not only to a proper 
extrapolation to vanishing volume fraction, but also to 
parameters characterizing the interaction of the parti- 
cles: Evidently, Eq. (jSJ refers only to equilibrium states. 
The analysis of I{q) obtained in the dilute regime may 
therefore yield to criteria that allow us to exclude non- 
equilibrium states. 

Arguments derived from first principles demonstrate 
that S(q) can be expanded into a series of powers in q"^ 

S{q) = S{0) + aq^ + (9) 

where the coefficient a is related to the second moment 
of the interaction potential V{r) and provides a measure 
for the range of interaction of the particles in suspension 
(see the discussion of this point in Ref. ^5^). In contrast 
to Eq. ||HJ), this relation holds for all possible volume frac- 
tions. Therefore Eq. Q provides another check for the 
data. In particular, it is highly useful to exclude scatter- 
ing intensities that have been affiicted by the presence of 
aggregates 0. 

C. Polydispersity 

As mentioned above, previous work has shown that 
Laponite particles are not monodisperse. The effect of 
polydispersity must therefore be taken into account and 
the intensity becomes an appropriate average over all 
species i 

f: c,M,P,{q) 

I{q) = K'-^^ SM{q) = KcMMQ)SM{q) , 

(10) 

where Pi (q) is the normalized form factor of a platelet of 
species i [Eq. ©], Sm is the "measured structure factor" , 
and s is the total number of species. Moreover, we have 
Ei=i '^i ~ ^- "The weight-average molecular weight 
of the platelets follows from this as 

M^ = ^4 • (11) 

E c» 



The normalized form factor of the multicomponent sys- 
tem Pwiq) can be expressed through the form factors 
Pi (g) of platelets of species i by 

Pu,{q) = —^ , (12) 

E 

i=l 

again with Pw{0) — 1. Obviously, all quantities discussed 
in the previous section as e.g. Bapp and deff must be 
replaced by appropriate mean values. 

D. Multicomponent interaction site model 

A quantitative understanding of correlations and inter- 
actions between various colloidal species can be achieved 
usin g th e well-established techniques of liquid-state the- 
ory |4a|. The PRISM ("polymer reference interac- 
tion site model") theory, originally designed for the 
study of polymer solutions and melts ^51j . has been re- 
cently extended to investigate various c harg ed and un- 
charged colloidal systems [ll El El EHIH El El , 
and mixtures of neutral spherical colloids and polymers 
60, 61,|6lE3l, as well as like- and oppositely-charged 
colloids and polyelectrolytes [3l0|. The key theoretical 
problem to be addressed here is that of an interpretation 
of experimental scattering data on very dilute Laponite 
clay suspensions within the framework of PRISM. The 
main difference of this work from the previous theoretical 
studies on monodisperse ^3 and bidisperse 's^ platelet 
suspensions is that we take into account polydispersity 
in the size of the platelets by using a multicomponent 
PRISM integral equation theory. We consider a multi- 
component system involving s species of charged platelets 
with number densities n;, where 1 < i < s. Each platelet 
of species i contains Ni equivalent interaction sites. The 
interaction potential between sites on particles of species 
i and j, carrying the charges ZiC and zje, will be of the 
generic form: 



Uij{r) ^ — — exp(-KDr), (13) 
er 

where Xd = is the usual Debye screening length. The 
system under investigation is an aqueous dispersions, but 
in view of the mesoscopic scale of the platelets, the sol- 
vent is modelled as a structureless dielectric continuum 
providing a macroscopic permittivity e. Any microscopic 
counterions or ions from added electrolyte will be consid- 
ered at the linear response (or Debye-Hiickel) level, i.e., 
they will screen the electrostatic potential due to the in- 
teraction sites on the charged platelets on a scale given 
by the Debye screening length. The underlying Born- 
Oppenheimer-like assumption entails that the charge dis- 
tribution on the mesoscopic particles does not contribute 
to screening. 
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FIG. 1: Absolute scattering intensities of Laponite suspen- 
sions normalized to their volume fraction (j> indicated in the 
graph. The data have been obtained by static light scattering 
in the region of smallest q-values (7.9 x 10~^ — 2.3 x 10~^ 
nm~^). Data at higher q were taken by using'the beamline 
ID2 (2.15 X 10"^ to 0.67 x 10"^ nm"^) and by a conventional 
Kratky-camera (0.07 to 6 nm~^). The solid line refers to the 
intensity extrapolated to vanishing concentration. 



PRISM is based on the assumption that all direct 
correlation functions between sites on pairs of different 
platelets are identical. This leaves a total of s(s + l)/2 
independent direct correlation functions Cy {q) . Corre- 
sponding total correlation functions hij (q) are defined by 
averaging over the NiNj correlation functions between all 
pairs of sites i and j on two platelets of the same or dif- 
ferent species. The two sets of correlation functions are 
related by the generalized Ornstein-Zernike equations of 
the PRISM, which in Fourier space read 



hij{q) = uji{q)cij{q)uji{q) + u!i{q)^Cii{q)niNihij{q) , 

1=1 

where NiUi is the number density of interaction sites of 
platelets of species i and uJi{q) = NiPi{q). The set of 
s(s -|- l)/2 independent Ornstein-Zernike relations must 
be supplemented by as many closure relations between 
each pair of total and direct correlation functions. Here 
we adopt the Laria-Wu-Chandler closure {6^] which has 
recently been used for charged polymers and colloids [l^ 



(15) 



where the asterics ★ denotes a convolution product. In 
addition, steric effects are taken into account using the 
Percus-Yevick approximation |46j |. 

The closed set of equations (|14|) and (|15|l are solved 
numerically by a standard iterative procedure to obtain 




FIG. 2; Extrapolation of the measured scattering intensities 
I{q) to vanishing concentration according to Eq. (0 and J^J. 
Parameter of the curves is the magnitude q of the scattering 
vector. The respective values are (in nm~^): quadrangles: 
0.0215; crosses: 0.037; open circles: 0.0526; open triangles: 
0.0681; filled circles: 0.0759; filled trangles: 0.0837. 



partial structure factors 

Sijiq) = UiNiUJiiq) + ninjNiNjhij{q) . 



(16) 



Scattering experiments measure the intensity of scattered 
radiation as a function of momentum transfer. To a mul- 
tiplicative constant this quantity is the same as the sum 

s 

of the partial structure factors ^ Sij (q) which is de- 

noted as the static structure factor. Most experimental 
papers on polyectrolytes, however, report the ratio of the 
static structure factor to the form factor. Therefore we 
compare in the following section the theoretical predic- 
tion for the function 



Siviiq) 



E S^J{q) 

s 

n,N,U!i{q) 



(17) 



to experimental data. In keeping with the experimental 
terminology we refer to S]\{ (q) as the measured structure 
factor. 



IV. RESULTS AND DISCUSSION 



Static scattering experiments: consistency 
check 



As mentioned in the Introduction, the present investi- 
gation focuses on the dilute regime in which no gelation 
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FIG. 3: Intensity Io{q)/(j> extrapolated to vanishing volume 
fraction <p. The data have been obtained by static light scat- 
tering in the region of smallest g- values (7.9x 10"'' — 2.3x lO"'^ 
nm~^). Data at higher q were taken by using the beamline 
ID2 (2.15 X 10"^ to 0.67 x 10"^ nm"^) and by a conventional 
Kratky-camera (0.07 to 6 nm~^). The solid line marks the fit 
according to Eq. 



occurs. However, no concentration is small enough to dis- 
regard the influence of mutual interaction of the platelets, 
i.e., of S{q). Moreover, the discussion in section lTlIBI has 
demonstrated that scattering data in this regime furnish 
valuable data as Bapp that may be used for a consistency 
check. As a first step in this analysis scattering data are 
therefore taken at small concentration (cj) < 0.16 %) and 
carefully extrapolated to vanishing concentration. Since 
we aim at the limits for vanishing scattering angles for 
all quantities under consideration, the g-range accessible 
by SAXS was not sufficient. Hence, static light scatter- 
ing was used to obtained data at smallest g-values pos- 
sible in order to supplement the analysis by synchrotron 
SAXS. For the highest scattering angles the data ob- 
tained by synchrotron SAXS were supplemented by in- 
tensities taken by a conventional Kratky-camera. This 
device can take data up to 6 nm~^. 

In all cases absolute scattering intensities have been 
obtained. For light scattering, the constant K in Eq. Q 
was found to K = (4.99 ± 1) 10"* cm^ Mol/g^ from the 
measured refractive increment [dfi/dc — 0.083 ± 0.008 
mL/g). In case of SAXS the respective optical constant 
K in Eq. Q was found to be if = 0.0445 Mol/g^. All 
absolute intensities agreed within the present limits of 
error (±5 %). 

Previous investigations by Bonn et al. and by 

Nicolai and Cocard [o^ have shown that even dilute 
aqueous suspensions of Laponite may contain aggregates 
that must be filtered off, in particular for light scattering 
experiments. Hence, all suspensions have been carefully 
filtrated prior to use. Moreover, special care was exerted 
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FIG. 4: Apparent second virial coefficient Bapp (cf. Eq. Q) 
of Laponite suspensions in the dilute regime, the filled circle 
denotes the value obtained by static light scattering whereas 
the other points have been obtained by SAXS. The effective 
diameter d^fj deriving from this fit results to 48 nm. 



in preparing the solutions by carefully stirring for a pro- 
longed time. If the time for the process of dissolution 
was chosen to be short or without vigorous stirring, ag- 
gregates resulted that showed up in the region of smallest 
angles. Evidently, the pH 10 and the ionic strength are 
further parameters of central importance. Most of the 
data to be discussed here were taken from solutions in 
which the ionic strength was set to 0.001 mM by NaCl. 

Figure ^ displays the data obtained at low ionic 
strength for various volume fractions (p calculated from 
the weight concentrations of the solute and its density 
in solution (2.42 ± 0.12 g/cm^). First of aU, the data 
obtained from SLS fit well together with the data ob- 
tained from SAXS. Moreover, the infiuence of the struc- 
ture factor S{q) cannot be disregarded even at the lowest 
volume fraction. This is due to the electrostatic stabiliza- 
tion that must lead to an effective diameter of interac- 
tion de// which is larger than the largest dimension of 
the platelets. Figure |21 demonstrates that these data can 
be safely extrapolated to vanishing concentration by use 
of Eq. (O and ijSll. Plots of c/I{q) vs. c lead to straight 
lines. The resulting intercepts are given as solid lines in 
Fig. n Even at the smallest concentration (</) — 0.02 %) 
the scattering intensities deviate from the value resulting 
from this extrapolation. Only at higher scattering angles 
the influence of mutual interaction is absent in full accord 
with Eq. ©. 

Figure 131 displays the intensity extrapolated to vanish- 
ing concentration for the entire g-range available from the 
combination of three devices, namely static fight scatter- 
ing, synchrotron radiation, and a conventional Kratky- 
camera. The rather flat region of small angles is due 
to the fact that Io{q) must follow Guinier's law (see 
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FIG. 5: Consistency check of the structure factor 5'i\f(g) in 
the region of small g- values according to Eq. ^ and 



Eq. ^) at small values of q. From a Guinier plot of 
these data (not shown) we obtained a weight average of 
the platelets Mu, = 930 ± 190 kg/Mol and a radius of gy- 
ration Rg = 13.4 nm. At intermediate q the slope of —2 
is seen predicted by Eq. ©. Finally, at highest 9- values 
the scattering intensity is greatly diminished because of 
the first minimum of P{q) which is located at g = 7.1 
nm~^. 

The form factor P{q) of a thin platelet can be eas- 
ily calculated according to Eq. ©. From this fit it 
becomes obvious that polydispersity must be included 
to give a satisfactory fit. This is due to the fact that 
Eq. © leads to a wavy from of Io{q) for monodis- 
perse platelets which is smoothed out by polydisper- 
sity. The solid line shows the resulting fit obtained for a 
Schulz-Zimm-distribution [btI 0| with a polydispersity 
expressed through M^/Mn = 1.5. The weight-average 
radius is given by 10.5 nm, the thickness is found to 
be 0.9 nm, and the calculated weight average [Eq. ((TT)l ] 
is in agreement with the measured value. We reiterate 
that Laponite platelets exhibit a appreciable polydisper- 
sity that can be deduced from this analysis. Obviously, 
this must be taken into account when modelling the in- 
teraction of the particles at finite volume fraction. 

The slope of the curves shown in Fig. El can now be 
evaluated to yield the apparent second virial coefficients 
Bapp defined through Eq. 0. Figure 0] displays the re- 
sulting data and a fit according to Eq. /.From this 
fit the effective diameter de// was obtained to 48 nm 
which is considerably higher than the largest dimension 
of the platelets (see above). This can be traced back to 
the electrostatic repulsion that is already operative be- 
fore the platelets touch each other. Hence, deff > 2R^ 
and Bam cannot be calculated from the dimensions of 
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FIG. 6: Structure factor Sm{<i) determined experimentally for 
various volume fractions (open circles: 0.02 %; open squares: 
0.04 %; open triangles: 0.08 %; solid circles: 0.12 %; solid 
squares: 0.16 %). The solid lines represent the correspond- 
ing structure factors computed according to Eq. H17|l and as- 
suming a purely repulsive interaction between the platelets 
[Eq. Cnj]- The inset shows the consistency check of the ex- 
perimental data according to Eq. @. 



the platelets in solution. The inset gives the same plot 
but with the abscissa scaled by (pdgff as suggested by 
Eq. ||SJ). A similar plot is displayed in Fig. |S1 showing 
the extrapolation to vanishing concentration as in Fig. |31 
Figure demonstrates that approximately parallel lines 
result from this scaling of the abscissa. Thus, the data 
obtained at low concentrations follow Eq. ©. This in- 
dicates that these solutions fully qualify as equilibrium 
states. 

The intensity referring to cj) — > may now be used 
to determine the structure factor S{q) for nonzero val- 
ues of (f). The foregoing discussion has clearly revealed, 
however, that the Laponite particles exhibit an apprecia- 
ble polydispersity. The structure factor obtained from 
this division should therefore be treated as a " measured" 
quantity S'a/('?) that contains explicitly a contribution 
of the size distribution. Figure [HI displays the result- 
ing SM{q) obtained for the low ionic strength of 0.001 
mM. The inset shows the consistency check according to 
Eq. It underscores again the fact that we are deal- 
ing with equilibrium states. Figure [HI shows that SMiq) 
is considerably lowered at vanishing scattering angle and 
exhibit a shallow maximum around q — 0.11 nm"-'^. As 
mentioned above, most of the investigations using scat- 
tering methods have been conducted near the concentra- 
tion of gelation, i.e., at much higher concentrations, and 
the structure factor has not been extracted. Hence, we 
cannot compare the present data directly with similar 
measurements from literature. 
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FIG. 7: The structure factor at zero scattering vector 
Sm (q = 0) (lower symbols) together with the structure factor 
at the main peak SM{q = q,nax) (upper symbols) for vari- 
ous suspensions: Laponite at five volume fractions: 0.02 %, 
0.04 %, 0.08 %, 0.12 %, 0.16 % from left to right (squares); 
polystyrenesulfonate of length L = 6.8 nm at two concentra- 
tions: 0.1 mol/L and 0.2 mol/L from left to right (up trian- 
gles 113]); polystyrenesulfonate of length L = 40 nm at two 
concentrations: 0.1 mol/L and 0.2 mol/L from left to right 
(down triangles 0); DNA of length L = 57 nm at 0.05 mol/L 
(circles "s^); DNA of length L = 380 nm at 0.05 mol/L (dia- 
monds 83]); proteoglycan at two different salt concentrations: 
0.05 mM and mM from left to right (plus-symbols fTTf ): to- 
bacco mosaic virus at five concentrations: 0.11 mg/mL, 0.27 
mg/mL, 0.43 mg/mL, 1.05 mg/mL, 2.07 mg/mL from left to 
right (crosses |84|]'l. 



The results for the structure factor as calculated ac- 
cording to Eqs. H13I) - H17|l are compared to the exper- 
imental scattering data in Fig. The adjustable pa- 
rameter related to the potential of interaction is the 
effective charge of the platelets which is much smaller 
than the structural charge of the platelets, due to coun- 
terion "condensation" 0, 0, [Z^- It is the strong re- 
duction of the bare charge which allows the use of lin- 
earized Poisson-Boltzmann theory, leading to the simple 
exponential Debye screening of the Coulomb interaction 
[Eq. (|13|l ]. In the calculations we use a fixed surface 
charge density which is defined as the charge of a platelet 
of species i divided by the area Ai = 7ri?| of the face 
of the platelet. The calculated structure factors agree 
with the experimental data for small scattering vectors, 
while pronounced deviations are visible in the intermedi- 
ate scattering vector regime. Particularly, the theoreti- 
cal calculations predict considerably smaller peaks of the 
structure factor than observed experimentally. Increas- 
ing the surface charge density, as compared to the one 
used in Fig. El leads to an increase of the height of the 



FIG. 8: Slow aggregation of Laponite particles at higher ionic 
strength (0.005 mM). The intensity Jo (?)/</!> was extrapolated 
to vanishing volume fraction (j) (cf. the discussion of Fig. . 
Note the much higher intensity obtained at small scattering 
vectors as compared to the case of the lower ionic strength 
(0.001 mM) in Fig. E] 



peaks of the calculated structure factors and to a decrease 
of the structure factors at small scattering vectors g — > 0. 
We emphasize that neither using different model parame- 
ters (surface charge density, Debye screening length, and 
size of the platelets) nor using a different distribution 
function (Gaussian, log-normal, and Tung distribution) 
lead to an agreement with the experimental data. 

Since structure factors and pair correlation functions of 
polyelectrolyte solutions calculated within the same theo- 
retical framework are in good agreement with experimen- 
tal data and computer simulations [s^ |^ JSj M^i 56], 
it is worthwhile to compare scattering data for poly- 
electrolytes with the data obtained for Laponite. Fig- 
ure [7| displays Sm{(1 — 0) which is proportional to the 
isothermal compressibility together with Sm{<1 = Imax) 
for various polyelectrolytes and Laponite, where qmax is 
the absolute value of the scattering vector at the main 
peak. The structure factors of the polyelectrolytes ex- 
hibit a main peak Sm{<1 = Qmax) > 1 and a small value 
Sniq = 0) < 0.11 for small scattering vectors, while the 
structure factors of Laponite are charcterized by a main 
peak and a rather large value 0.7 < Sm{(1 = 0) < 0.95 
for small scattering vectors. The enhancement of the 
small angle {q 0) value of SM{q) of Laponite as com- 
pared to the polyelectrolytes signals increased density 
fluctuations. Since the polyelectrolyte solutions remain 
liquid-like even at higher concentrations, the observed 
qualitative different behavior of the structure factors of 
Laponite may be considered as indicative of the sol-gel 
transition at higher concentrations. Notice that the val- 



ues of Sufiq = 0) for the Laponite suspensions under 
considerations are smaller than 1 and decrease with in- 
creasing concentration in contrast to strong small q up- 
turns which have been observed experimentally for clay 
suspensions in the gel phase ,18..„ .20.. ..2L ■7L,72> .73,, ,74 .75| . 
polyelectrolyte gels and mixtures [23|) polystyrene- 
sulfonate ion exchange resins , and low ionic strength 
polyelectrolyte solutions J9_, 80, 81). These strong up- 
turns signal strong concentration fluctuations indicative 
of aggregation, or spinodal instability reminiscent of the 
behaviour observed in recent c omp uter simulations [s^ 
and PRSIM integral equations [Sg of mixtures of oppo- 
sitely charged particles. 

A final point to be discussed at this point is the stabil- 
ity of Laponite suspensions at higher ionic strength. Fig- 
ure |H1 displays the intensity obtained for an ionic strength 
of 0.005 mM. The same procedures have been applied as 
in case of the data discussed above. The intensities have 
been extrapolated to vanishing concentration in order to 
remove all alterations due to S]\i{q). The dashed line 
marks the form factor obtained previously for 0.001 mM 
as described above. It is obvious that the intensity mea- 
sured at small scattering angles is considerably higher 
than in case of the lower ionic strength. We interpret 
this as the onset of a slow coagulation process. Hence, a 
salt concentration as low as 0.005 mM already suffices to 
lead to the built-up of aggregates visible in the scattering 
experiment at small q. This points clearly to the limited 
stability of the Laponite platelets in solutions which will 
be further discussed below. 



B. Effective interaction potential 

The results of the preceding subsection indicate that 
the disagreement between the experimental and theoret- 
ical results may be due the purely repulsive screened 
Coulomb potential used in the calculations [Eq. H13|) ]. 
Recently Ruzicka et al. 0| have suggested that at low 
concentrations the effective interaction between Laponite 
platelets are characterized by a competition of long-range 
Coulomb repulsion and short-range van der Waals attrac- 
tion similar to colloidal systems and protein solutions 
HE Ell' Moreover, the charge density is higher on 
the face than on the rim of a platelet, which leads to a 
modification of the interaction potential between sites on 
particles i and j [Eq. Ijl^fl ] within the computationally 
very demanding multicomponent PRISM theory. Here 
we attempt a mesoscopic coarse graining, whereby par- 
ticles act via an effective potential V{r). To this end we 
solve numerically the Ornstcin-Zernike equation 

h{q) ^ c{q) + nc{q)h{q) (18) 
together with hypernetted-chain closure relation 

h{r) = \n[h{r) + 1] + c{r) + {kBT)-^V{r) . (19) 




0.1 0.2 0.3 
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FIG. 9: Structure factor Sm^q) determined experimentally for 
various volume fractions (open circles: 0.02 %; open squares: 
0.04 %; open triangles: 0.08 %; solid circles: 0.12 %; solid 
squares: 0.16 %). The solid lines represent the results of the 
integral equation theory [Eqs. 1)18^ - l|2()|l ] using the effective 
potential V{r) shown in the inset. 

These equations follow from Eqs. (|14|l and H15|) for s = 1, 
LOi{q) = 1, and Ni = \. The calculated structure factor 

S{q)^\ + nh{q) (20) 

is compared in Fig. ^ with the experimental data. For 
all volume fraction, good agreement with experiment is 
achieved with a density-independent effective potential 
V{r) that is attractive for small distances but repulsive 
for larger distances (see the inset in Fig. 0). We em- 
phasize that the effective potential fulfills the stability 
condition drr^V{r) > discussed by Ruelle 89], 
and provides a representation of the underlying many- 
body interactions in the system. Thermodynamic prop- 
erties can be calculated from the compressibility relation 
S{q = 0) = nksTnT, where kt is the isothermal com- 
pressibility |4(t| |. The osmotic pressure is evaluated as 
P = UbT J dnS{q ~ 0)^^. Moreover, the correlation 
functions obtained from the Ornstein-Zernike equation 
can be used as input for a Langevin equation for the 
dynamic structure factor [5^1 allowing one to study dy- 
namical properties. 

Finally, it is instructive to rationalize our results us- 
ing simple arguments. The repulsion due to the over- 
lap of the counterion diffusive layers can be replaced 
by a renormalized hard-core repulsion with an effec- 
tive diameter equal to 2(i?u, + \o) = 41 nm, where 
the Debye screening length has been approximated by 
\b = (STrZans)"-"^/^, where Ib = 0.714 nm is the Bjer- 
rum length and 2ns is the particle number density of 
salt ions at ionic strength 0.001 mM. This effective di- 
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ameter is somewhat smaller than de// = 48 nm ob- 
tained from Eq. (jS)). Another length of interest is the 
estimated average distance between two platelets defined 
as dh = where = ttR^^L is the average 

particle volume. For cj) — 0.02 % we find dh = 116 nm 
and for (j) = 0.16 %, dh ~ 58 nm. As is apparent from 
Fig. O the correlation length associated with the main 
peak of the structure factor, dc = 'ii^/qmax ~ 53 nm 
is independent of the volume fraction in contrast to re- 
sults for Laponite at very low ionic strength (0.0001 M; 
[T^I)- The correlation length dc for the samples shown 
in Fig. 1^1 is related to the effective diameter def / rather 
than to average distance between two platelets defined 
as dh- However, the fact that deff is comparable to 
dh for (j) — 0.16 % indicates that the sol-gel transition 
occurs already at a volume fraction slightly larger than 
4> — 0.16 % for which de// ~ dh, in agreement with our 
experimental findings. The main peak at the scattering 
vector Qmax for gel systems reflects the short interparticle 
distance between connected particles of effective diame- 
ter. The pronounced excess scattering at small values 
of q is due to correlations associated with the network 
structure. 



nitude q of the scattering vector so that reliable data 
of the structure factor SMiq) could be obtained. Spe- 
cial care was taken to show that the present data re- 
fer to true equilibrium states. The form factor obtained 
from the scattering intensity extrapolated to vanishing 
volume fraction agrees with the calculated one for poly- 
disperse platelets (see Figure ISJ- The structure fac- 
tor SM{q) is characterized by a main Coulomb peak 
Sniq — qmax) > 1 and a value 0.7 < S'm(9 = 0) < 0.95 
for small scattering vectors corresponding to a rather 
large isothermal compressibility, in contrast to various 
polyelectolyte solutions (see Figure [Tj) . The analysis of 
Sm {q) which includes the marked effect of polydispersity 
demonstrates unambiguously that the Laponite platelets 
interact via a potential that is attractive on short dis- 
tances but repulsive on longer distances (see inset of Fig- 
ure 0). This clearly points to the role of short-range at- 
traction for the process of coagulation at higher volume 
fractions. Moreover, it explains that higher concentration 
of added salt lead to aggregation because of the increased 
screening of the repulsive Coulomb interaction between 
the Laponite platelets (see Figure [HJ. 



V. CONCLUSION 

A precise analysis of the static scattering intensity of 
dilute aqueous Laponite solutions has been presented. 
The combination of static light scattering with small- 
angle x-ray scattering leads to a wide range of the mag- 
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